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ABSTRACT 

Eukaryotic translation termination results from the 
complex functional interplay between two release 
factors, eRF1 and eRF3, in which GTP hydrolysis 
by eRF3 couples codon recognition with peptidyl- 
tRNA hydrolysis by eRF1. Here, we present a cryo- 
electron microscopy structure of pre-termination 
complexes associated with eRF1»eRF3»GDPNP at 
9.7 -A resolution, which corresponds to the initial 
pre-GTP hydrolysis stage of factor attachment and 
stop codon recognition. It reveals the ribosomal 
positions of eRFs and provides insights into the 
mechanisms of stop codon recognition and trigger- 
ing of eRFS's GTPase activity. 

INTRODUCTION 

Termination occurs when a stop codon enters the riboso- 
mal A site and consists of stop codon recognition followed 
by peptide release, which involves the nucleophilic attack 
of a water molecule on the P-site peptidyl-tRNA in the 
ribosomal peptidyl transferase center (PTC). Eukaryotic 
termination is mediated by two directly interacting release 
factors: eRFl, which is responsible for stop codon recog- 
nition and triggering peptide release, and eRF3, a GTPase 
that strongly stimulates peptide release by eRFl in a GTP- 
dependent manner (1-3). eRFl, in turn, stabilizes binding 
of GTP to eRF3 so that they form a stable ternary 
complex (4,5), and is required for eRF3's ribosome- 
dependent GTPase activity (6). eRFl has omnipotent 
decoding capacity and recognizes all three stop codons. 



It comprises N-terminal (N), middle (M) and 
C-terminal (C) domains (7). The rigid core of domain 
C contains a flexible insertion forming a mini-domain 
(8). Domain N is involved in stop codon recognition. 
Although the mechanism by which eRFl responds to all 
three stop codons is not clear, extensive mutational and 
genetic analyses identified the essential role in this process 
of GTS3i_33 (human numbering), TASNIKS58-64 and 
YxCxxxF 125-1 31 motifs located at the apex of the N- 
domain (7,9-18). eRFl's domain M contains the universal 
GGQ loop, whose placement into the PTC causes re- 
arrangement of rRNA, allowing a water molecule to 
enter and induce peptide release (1,7,19-22). 

eRF3 consists of the essential C-terminal region 
comprising GTP-binding (G) domain and (3-barrel 
domains 2 and 3, which are homologous to elongation 
factors EF-Tu and eEFlA (23), and a non-conserved N- 
terminal region that is not essential for termination (24). 
eRFl's C and M domains interact with eRF3 (15,25). 
eRFl and eRF3 bind to the pre-termination complexes 
(pre-TCs) as an eRFl»eRF3»GTP ternary complex, but 
peptide release does not occur until eRF3 hydrolyzes 
GTP. GTP hydrolysis releases eRFl's domain M from 
eRF3, which enables its GGQ loop to enter the PTC 
and trigger peptide release. 

Although the structures of individual eRFs have been 
determined (7,15,23), the structural basis for key steps 
in termination, such as stop codon recognition or trigger- 
ing of eRF3's GTPase activity, remains unresolved. Our 
previous cryo-EM structure of the eRFl»eRF3»GDPNP- 
bound pre-TC had a resolution of 18 A (26), which did not 
allow accurate modeling of the ribosome-bound factors 
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and therefore could not address such questions. The 
advent of a more powerful classification algorithm, 
RELION (27), prompted us to re-examine the original 
data set (26), resulting in a much higher resolution 
(9.7 A), which enabled us to propose a more detailed 
model of the interactions pivotal for peptide release. 

MATERIALS AND METHODS 

Data processing 

Mammalian pre-TCs were assembled in vitro on a deriva- 
tive of (3-globin messenger RNA (mRNA) encoding 
MVHL tetrapeptide followed by a UAA stop codon, 
using purified rabbit 40S and 60S ribosomal subunits, ini- 
tiation (2,3,1,1A,4A,4B,4G,5 and 5B) and elongation (lA 
and 2) factors and aminoacylated tRNAs (26). 
eRFl»eRF3»GDPNP-bound pre-TCs were formed with 
full-length eRFl and eRF3 lacking the non-essential 
N- terminal 138aa (26). A cryo-EM data set of the 
eRFl»eRF3»GDPNP-bound pre-TCs comprising 
195 432 particles (26) was aligned and refined using the 
RELION autorefine procedure (27). It resulted in a 
9.1 -A structure [gold standard Fourier shell correlation 
(FSC) = 0.143] showing fragmented densities in the 
intersubunit space, indicative of heterogeneity. The data 
set was then subjected to RELION classification (28) 
starting with 10 classes (k = 10). Of the 10 classes, class 
8 (Figure 1) was well populated (48 973 particles, 25.1%), 
with well-defined additional masses of density in the 
intersubunit space corresponding to eRFl, eRF3 and P- 
site tRNA as previously described (26). The particles from 
this class were isolated, and the structure was first refined 
with RELION autorefine procedure. This yielded a struc- 
ture at 9.7 A (gold standard FSC = 0.143), but the density 
exhibited strong stretching artifacts when filtered at the 
measured resolution, probably due to preferred orienta- 
tions of the particles on the cryo-EM grid. The aUgnment 
parameters were further refined using RELION classifica- 
tion with k = 1 and a T factor of four for three iterations 
at an angular spacing of 1.8°, followed by three additional 
iterations at 0.9°. The resulting reconstruction displayed 
the same resolution of 9.7 A (FSC = 0.143), but the 
quality of the structure was very much improved and the 
stretching artifacts were no longer visible. The statistics of 
aUgnment precision estimated by RELION were also sig- 
nificantly improved [improvement in terms of overall 
accuracy of rotations, 41%; overall accuracy of transla- 
tions, 43% and average Pmax at 1.8° angular step, 180% 
(Supplementary Table SI)]. 

Modeling of eRFl, eRF3, mRNA and P-site tRNA^^" 

The cryo-EM map segmentation was performed using the 
Segger (29) plug-in in UCSF Chimera (30). Models of 
human eRFl and eRF3 were built based on several experi- 
mental structures. Atomic coordinates of eRFl and eRF3 
domains 2 and 3 were taken from their crystal structure 
(15) (PDBID: 3E1Y), and the C-terminal mini-domain 
structure of eRFl (residues 326-373) was modeled based 
on its nuclear magnetic resonance (NMR) structure (8) 
(PDBID: 2KTU). Most of eRF3 G-domain (residues 



317-439, 207-216) was modeled by homology based on 
the crystal structure of Saccharomyces pombe eRF3»GDP 
(23) (PDBID: 1R5B) because of its significantly higher 
resolution over the S. pombe eRF3»GMPPNP crystal 
structure (23) (2.35 versus 3.20 A, respectively). Only the 
GTP-binding pocket of the G-domain (residues 217-231) 
was modeled based on the crystal structure of 
eRF3#GMPPNP (23) (PDBID: 1R50). Missing residues 
and the switch region (residues 232-316) in eRF3 G- 
domain were modeled by homology based on the archeal 
elongation factor la (aEFla) taken from the eRFl»GTP- 
bound aEFla crystal structure (31) (PDBID: 3VMF). The 
arrangement between eRF3's domains II and III and G- 
domain was derived from the latter crystal structure as 
well. The P-site tRNA^^^ was modeled from its crystal 
structure in complex with leucyl-tRNA synthetase from 
Pyrococcus horikoshii (32) (PDBID: 1WZ2). The mRNA 
was modeled based on the mRNA of the crystal structure 
of the Thermus thermophilus 70S ribosome bound with the 
Q253P mutant form of release factor 2 (22) (PDBID: 
4KFK). All homology modeUng was done using SWISS- 
MODEL (33,34). 

Molecular dynamics flexible fitting 

The derived model was first rigid body-fitted into its cor- 
responding segmented cryo-EM density map using UCSF 
Chimera (30). Then, using the molecular graphics software 
VMD (35), the system was prepared for molecular 
dynamics flexible fitting (MDFF) (36) in expHcit solvent, 
a procedure which applies the cryo-EM map as an add- 
itional potential to the system — in this study the 
eRFl»eRF3»tRNA»mRNA cryo-EM segmented map — 
thus comprising only the molecules to be simulated: 
eRFl, eRF3, mRNA and P-site tRNA^'^. The MDFF 
system was embedded in a box of TIP3P water molecules 
with an extra 12- A padding in each direction. The system 
was neutralized by potassium ions, and an excess of 
~0.2M KCl was added. The simulated system was 
prepared using CHARMM force field parameters 
[Combined CHARMM All-hydrogen topology file for 
CHARMM22 proteins and CHARMM27 lipids (37)]. 
The system was energy-minimized in 600 steps in the 
molecular dynamics simulation package NAMD (38). 
After minimization, the fitting trajectories were run for 
400 ps, once the root-mean-square deviation and cross- 
correlation coefficient had stabilized. 

RESULTS AND DISCUSSION 

Structure determination of the mammalian 
eRFl#eRF3#GDPNP-bound pre-TC 

Classification of the cryo-EM data set comprising 195 432 
particles of the mammalian eRFl»eRF3»GDPNP-bound 
pre-TC formed on mRNA encoding MVHL tetrapeptide, 
followed by a UAA stop codon, yielded 10 classes with 
different factor and tRNA occupancies (Figure 1). One of 
these classes (25% of the particles) had P-site tRNA along 
with well-defined eRFl and eRF3 densities. The charac- 
teristic long variable loop of tRNA^^^ was apparent in the 
P-site tRNA, indicating that it was correctly programmed 
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with an A-site stop codon. The eRFl»eRF3»GDPNP- 

bound pre-TC class yielded a 9. 7- A reconstruction 
(Figure 2a), which allowed precise docking of each 
factor's domains. eRFl and eRF3 were modeled 
(Figures 2b-d) based on the crystal structures of eRFl, 
eRF3, eRFl-eRF32.3 and aRFl-aEFla (7,15,23,31) and 
the NMR structure of eRFl's C-domain (8). The model 
was refined using MDFF (36) with secondary structure 
constraints. 

Overview of the complex 

The position and conformation of ribosome-bound eRFl/ 
eRF3 are similar to those of aa-tRNA/EF-Tu (39) and the 
eRFl/eRF3 paralogs Dom34/Hbsl (40) involved in the 
dissociation of stalled elongation complexes (41,42), with 
eRFl's domain N binding to the decoding center, and the 
rigid core and mini-domain of eRFl's domain C forming a 
bridge between the P stalk and the beak of the 40S 
subunit. eRFl interacts with eRF3's domain 3 via its 
domain C, whereas eRFl's extended domain M contain- 
ing the GGQ motif is tucked between eRF3's domains 2 
and G. eRF3 is bound to the universal GTPase-associated 
center (GAC) of the ribosome, between the sarcin-ricin 
loop (SRL) on the 60S subunit, and helices (h) 5 and 14 
of 18S rRNA on the 40S subunit. No additional conform- 
ational changes were observed on the ribosome at the 
improved resolution, compared with those reported previ- 
ously (26): it is in a non-rotated state, the LI stalk in the 
open position (Figure 1) and the P stalk base shifted 
inward. However, what was previously interpreted as the 
result of a ribosomal rearrangement at the mRNA 
entrance (26) now appears more likely to be an mRNA 
bundle (Figure 2 and Supplementary Figure SI). 

Interaction between eRFl-N and the stop codon 

eRFl's N-domain reaches deep into the decoding center 
(Figure 2), establishing multiple contacts with the 40S 
subunit, including hi 8, h30, h31, h34 and h44 of 18S 
rRNA and ribosomal proteins rpS30e and rpS31e. R65/ 
R68 in the long a3 helix are in proximity to nucleotides 
1330-1331 of h34, consistent with these residues' strong 
influence on eRFl's ribosomal binding (11). In our struc- 
ture, mRNA occupies its normal position in the decoding 
center, and eRFl's N-domain interacts directly with the 
stop codon (Figure 3), primarily via the TASNIKS58_64 
motif, consistent with ultraviolet cross-linking of the 
stop codon to this element (10). The positions of the 
mRNA bases and amino acid side-chains cannot be 
determined precisely at the present resolution. However, 
the observed proximity of particular residues and bases in 
the fitted structures makes it likely that they interact with 
one another. The TASN residues appear to be close to the 
first position of the stop codon, I to the second and KS to 
the third. The proximity of K to the nucleotide immedi- 
ately downstream of the stop codon could potentially con- 
tribute to the mechanism by which this base influences 
termination efficiency (43,44). In addition to TASNIKS, 
three residues from GTS3i_33 and YxCxxxF 125-1 31 motifs 
are also positioned within interacting distance of the stop 
codon: G31 and T32 are close to the third and second 




Switch 1 

Figure 2. Overview of the complex, (a) Overview of the map showing 
the 40S (yellow), 60S (blue), eRFl (Magenta), eRF3 (red), P-site tRNA 
(green) and mRNA (coral), (b) Atomic model fitted into its density 
viewed from the 60S side with the 60S density removed, (c) Atomic 
model fitted into its density viewed from the 40S side with the 40S 
density removed, (d) Atomic model in the same orientation as in c, 
showing eRFl, eRF3, the P-site tRNA and mRNA path. 

base, respectively, and CI 27 to the bases in positions 1 
and 2. The positions of stop codon nucleotides relative to 
these three residues are consistent with the experimentally 
determined influence of T32 and CI 27 on the specificity of 
recognition of the nucleotide in the second position of the 
stop codon (14,15,17). Y125 and F131 are farther away, 
and their influence on termination (9,12,15,17) may there- 
fore reflect roles in stabiHzing the domain structure. 
Additional contacts with the pre-TC may be important 
for proper placement of the N-domain into the decoding 
site. Thus, the P-sheet supporting the GTS and YCF motifs 
interacts with H69 of 28 S rRNA, whereas the continuous 
density between helix a2 and the P-site tRNA (which is 
particularly strong between S46 at the tip of a2 and nucleo- 
tide 30 of the tRNA^^^ anticodon stem-loop) indicates 
likely contact between eRFl and P-site tRNA. 
Functional interaction between them has been proposed, 
and preferential binding to specific tRNAs could 
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Figure 3. Close-up of the eRFl N-domain and mRNA showing the TASNIKS amino-acids (blue), GTS loop (red) and YxCxxxF motif (green), and 
their positions relative to the approximate location of the stop codon bases represented as slabs. Codon positions are indicated in orange: first 
position, yellow: second position and blue: third position. 



contribute to the observed bias for certain codons 
preceding stop codons (45,46). Interestingly, the character- 
istic long variable loop of the P-site tRNA^^^ is also in 
proximity to the highly conserved P-site loop of rpLll on 
the 60S subunit and to h30 and rpS18 of the 40S subunit 
(Supplementary Figure S2). The P-site loop of rpLl 1 (L5 in 
bacteria) contacts the T-loop of tRNA in the P/P, as well as 
in neighboring intermediate states, and was therefore 
proposed to escort tRNA through the P site during trans- 
location (47-50). Notably, residues in the P-site loop, which 
are seen in proximity to the variable loop of tRNA in the 
present reconstruction, were shown to play a key role in 
binding of the P-site tRNA (51). The potential strengthen- 
ing of the contact with the P-site loop of LI 1 and additional 
specific contacts with h30 and rpS18 may contribute to the 
known higher affinity of tRNA^^"" to the P/P state, ac- 
counting for the increased stabiHty of ribosomal complexes 
containing P-site tRNA^^^ after peptide release (52). 

Additional stabilization of eRFl's interaction with the 
40S subunit appears to be provided by its mini-domain, 
which protrudes toward the beak where it interacts via its 
flexible loop with h33/h34 of 18S rRNA, or possibly with 
the N-terminal tail of rpS31e (Figure 4). 

Interaction between eRFl and eRF3 

The rigid core of eRFl's C-domain interacts with eRF3's 
domain 3 as observed in the eRFl-eRF32-3 crystal struc- 
ture (15), and with H44 of the P stalk (Figures 2 and 4). 
rpL12 is in proximity to eRFl's C-domain and domain 3 
of eRF3 and may contact them directly, but its scattered 
density suggests that they do not interact stably in this 
configuration. The inward rotation of the stalk base (26) 
is similar to that observed in Dom34/Hbsl -associated 
ribosomes (40). 

The topology of eRFl's M domain is now clearly vis- 
ible, folded back onto eRF3, with the GGQ loop close to 
the switch I region, similarly to the interaction between 
Dom34's central domain with Hbsl (40) and aRFl 



□ 




Figure 4. Close-up of eRFl-C and the mini-domain, showing the inter- 
action of eRFl-C with H44 of the P stalk, and the interaction of the 
mini-domain with the 40S beak. rpL12 and rpLPO not shown. 



domain B with aEFla (31) (Supplementary Figure S3), 
whereas a-helix 8 interacts with eRF3's domain 3, close 
to the conserved GRFTLRD613-619 motif (23) (Figure 5a). 

Structure of eRF3 bound to the ribosome 

eRF3's domains are well resolved: domains 2 and 3 
have the same orientation relative to eRFl as in the 



3414 Nucleic Acids Research, 2014, Vol. 42, No. 5 




Figure 5. (a) Close-up of eRF3 showing the eRFl-eRF3 model fitted in the electron density together with the human SOS ribosome model (PDBID: 
3J3A/B/D/F) [18S rRNA (purple), small subunit ribosomal proteins (beige), 28S rRNA (green), large subunit ribosomal proteins (hght blue)], (b) The 
70S ribosome/EF-Tu/tRNA/GDPNP from (39) compared with the SOS ribosome-eRFl-eRF3 in (a). 



eRFl-eRF32-3 complex (15), whereas the G-domain is 
positioned similarly to that of aEFla in the aRFl- 
aEFla complex (31) and to G-domains in other 
GTPases, such as EF-Tu (39), and contacts the GAC in 
a similar manner (Figure 5). The switch I region of eRF3's 
G-domain is well ordered and is tucked between the SRL, 
hl4 of 18S rRNA and the GGQ loop of eRFl-M. The 
SRL nt 4599-4600 and 4601-4603 are in proximity to 
DK273-274 in the switch I and GE326-327 in helix a5, re- 
spectively (Figure 5a). 

GTPase-activation mechanism 

There are numerous similarities in the functioning of 
eRFl»eRF3»GTP and aa-tRNA»EF-Tu»GTP complexes 
that suggest likely conservation of the GTPase-activation 
mechanism: hydrolysis of EF-Tu/eRF3-bound GTP is 



triggered by binding to the A site of aa-tRNA or of 
eRFl, respectively, followed by dissociation of the factor 
and accommodation of either the aa-tRNA acceptor end 
or the apical loop of eRFl's M domain in the PTC. In 
the case of EF-Tu, binding of aa-tRNA»EF-Tu»GTP to 
the ribosome induces a shift in domain 2, resulting in the 
interaction of a highly conserved (3-turn with h5 and 
disruption of the interaction of the 3 '-end of tRNA with 
the switch I loop, which leads to GTPase activation in a 
process that also involves direct participation of the SRL 
(39,53). A similar contact is established between h5 on the 
40 S subunit and the conserved KDxGT449_453 P-turn in 
eRF3's domain 2 (Figure 5a), and a C3302U substitution 
in the Saccharomyces cerevisiae SRL (equivalent to human 
C4603U) caused a termination defect (54). One notable 
difference is that there is additional density between hl4 
and eRFl's M domain (Figure 5a), which involves the 
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Figure 6. (a) Superimposition onto domain C of eRFl in the crystal structure (7) (PDBID: 1DT9) and bound to the pre-TC ribosome. The mini- 
domain is not displayed for clarity, (b) Comparison between the positions of eRFl-M in the eRFl-eRF3-ribosome structure, the eRFl-eRF32-3 (15) 
(PDBID: 3E1Y) and the eRFl (7) (PDBID: 1DT9) crystal structures with their domain C superimposed. Position of the P-site tRNA in the context 
of the pre-termination complex is shown to compare the distance necessary to be travelled for the GGQ loop to reach the peptidyl transfer center 
(PTC) with the amphtude of movement between the different conformations of eRFl observed. Only the N domain of eRFl in the pre-TC cryo-EM 
structure is shown for clarity. 
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conserved R192 that is important for GTP hydrolysis, but 
not for GTP binding or the eRFl/eRF3 interaction (15). 
A hypothetical mechanism of GTPase activation based on 
the presumed similarity with EF-Tu would be that off the 
ribosome, eRFl's M domain is too tightly packed onto the 
switch I region to allow GTP hydrolysis, but after 
eRFl»eRF3»GTP binds to the ribosome, the R192/hl4 
and (3-turn/h5 interactions together pull the M-domain 
from eRF3, loosening the switch I region. The P stalk 
may act as a gatekeeper in this regard, impairing the 
eRFl-M/hl4 and eRF3/h5 interactions until eRFl's 
binding to a bona fide stop codon is strong enough to pull 
the P stalk inward. Although eRFl can stimulate ribosome- 
dependent GTPase activity in the absence of an A-site stop 
codon (6), it is likely that GTP hydrolysis would be 
accelerated by the interaction between eRFl and a stop 
codon. 

In eRFl#eRF3#GDPNP-bound pre-TCs, the GGQ 
loop of eRFl is positioned far away from the PTC, and 
its accommodation in the PTC would require substantial 
rearrangement of eRFl's domain M. Superimposing the 
structure of eRFl in the eRFl»eRF3»GDPNP-bound 
pre-TC onto the crystal structures of individual or 
eRF32-3-bound eRFl (7,15) shows that domain M can 
undergo large hinge movements (Figure 6). If eRF3- 
bound eRFl adopts a strained conformation on binding 
to the ribosome similar to that of the A/T state of EF-Tu- 
bound tRNA (55,56), then dissociation of eRF3 following 
GTP hydrolysis would initiate spontaneous relaxation of 
eRFl resulting in accommodation of the GGQ loop in the 
PTC. Full entry of the GGQ loop into the PTC would 
then likely be facilitated by the high mobility of the 
apical GGQ loop in the thermal environment (57), the 
positive charge of this region promoting its retention 
in the PTC following entry. 



CONCLUSION 

Our structure revealed the topology of the interaction 
between eRFl's N-domain and the UAA stop codon, as 
well as the intricate network of interactions between 
eRF3's switch I region, eRFl's GGQ loop and the 
GAC. It suggests that GTP hydrolysis may be triggered 
once eRFl is properly bound to the STOP codon in a 
manner similar to that of EF-Tu. Now, atomic resolution 
structures of this complex bound to the three different 
stop codons may be required to achieve complete under- 
standing of how a single factor can recognize them all. 
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